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DISCLAIMER 



Although each program has been tested by its contributor, 
no warranty, express or implied, is made by the contributor 
•or any User's Group, as to the accuracy and functioning of 
the program and related program material, nor shall the fact 
of distribution constitute any such warranty, and no responsibility 
is assumed by the contributor or any User's Group, in connection 
therewith. 



Modifications or revisions to this program, as they occur, 
will be announced in the appropriate Catalog of Programs 
for IBM Data Processing Systems. When such an announce- 
ment occurs, users should order a complete new program 
from the Program Information Department. 
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COMMON USERS GROUP PROGRAM REVIEW AND EVALUATION 
© (fill out in typewriter, ink or pencil) 

Program No. ' ~ Date : 

Program Name: ; 

1. Does the abstract adequately describe what the program is and what Yes No_ 

it does ? 

Comment 

2. Does the program do what the abstract says ? Yes No_ 

Comment 

3. Is the description clear, understandable, and adequate? Yes No_ 

Comment 

4. Are the Operating Instructions understandable and in sufficient detail? Ye s No_ 
Commen t . 

Are the Sense Switch options adequately described (if applicable) ? Yes No_ 

Are the mnemonic labels identified or sufficiently understandable ? Yes No_ 

Comment 



5. Does the source program compile satisfactorily (if applicable) ? Yes No_ 

Comment 

6. Does the object program run satisfactorily?' Yes No_ 

C omm ent 

7. Number of test cases run . Are any restrictions as to data, 

size, 'range, etc. covered adequately in description? Yes No_ 

Comment v 

8. Does the Program meet the minimal standards of COMMON? Yes No_ 

Comment - • . 

9. Were all necessary parts of the program received? Yes No_ 

Comment 

10. Please list on the back any suggestions to improve the usefulness of the program. 
These will be passed onto the author for his consideration. 



Please return to: Your Name 

Company 

Address 



IBM Corporation Company 
Program Information Department 
40 Saw Mill River Road 



Hawthorne, New York 10532 Users GrounTvxle" 

Attn: PREP FORM COORDINATOR . . Users Group Code 



THIS REVIEW FORM IS PART OF THE COMMON ORGANIZATION'S PROGRAM REVIEW AND 
EVALUATION PROCEDURE. NONME MBERS ARE CORDIALLY INVITED TO PARTICIPATE 
IN THIS EVALUATION. 
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PROGRAM DESCRIPTION 



This prog/cam. is as nearly as possible a general purpose program for the 
frequently encountered curves y - b*c/(x-a) , y « ab*, and y • ax* 1 . It in 
really three subprograms, and manual branching Tia 1620 console is neces- 
sary to ■ore from one subprogram to another. All output is Tia the con- 
sole typewriter. All input is Tia the 1622 card reader or console type- 
writer (OPTION 5 only). For each sub pr ogram the options listed below are 
available to the user. 

0PTI0M 1 

Fit by Least-Squares Method a corre of desired form to a set of empirical 
data. When Op 1 is chosen, batch processing is possible. Normally batch 
p roc es sing would be used when experience dictated fona of curre to be used 
for several sets of data. 

OPTIOH 2, "2a" 

Op 1 pins standard error and plot back of as many points as desired. In- 
terpolated values furnished also. Batch processing is not possible, be- 
cause two passes through the card re ad er mast be made with the sane set of 
observed data. A third pass is used for interpolation. 

OPTION 3 

. Mo interpolation; otherwise sane as Op 2. 

OPTIOH 4 

Op 1 pins plot back of a few select points and/or interpolation. Ho stan- 
dard error computation. This is a two-pass option with a t h r e e pass func- 
tion, and careful reading of operating instructions is recommended. 

OPTIOH 5 

Function evaluation by keying-in parameters at console typewriter. Use of 
this feature permits evaluation of any function of specified form, provided 
arithmetic and function t alnes n— In within bounds of FORTRAN floating 
subroutines. Fea tu re especially useful when interpolation is "after- 
thought". 

METHOD OF COMPUTATION 
LKAST 9JMBES FIT: (HYPERBOLA) 

y - b+o/(x-e)— ^(y-b)(x-a) - c -^xy-bx-ay+ab - c =^xy - (c-ab)+bx+ay. 
Least squares method is used to fit last f on to set of n ob ser ved data 
(xo,y<>)- Three normal equations in the three unknowns, Tc-ab) , b, and a, 
result, lie solve these simultaneously: 

nCc-ab^b^x,, + a^y^ -X*o7o W 
o-l o»l o»l 



1 



2 



(c-ab)^xo + bUxo 2 + &Hx^j =Xx 2 y (2) 

(c-ab)2Ty + bUxoyo + aZIyo 2 =23oy 2 (3) 

Third equation was solved for (c-ab), and result was substituted in (1 ) and 
(2) from which a was then eliminated. Resulting expression for b required 
two FORTRAN statements. Value a obtsined by formula resulting from substi- 
tution of b in modified second equation. Value (c-ab) obtained by substi- 
tution of a and b into formula given by initial, solution of (3) for (c-ab). 
Subsequent solution for c_ is imr.ediate. 



LEAST SQUARES FIT: (EXPONENTIAL) 



y = at^^ln y = In a + (In b)x. Let In y = y', In a = A, In b = E. Hence 

y' = A + Bx. To fit line cf this form to a set of n observed data (xo,y i 

we perform the transformation y c ' = In y G . Then we solve simultaneously 
for A and B the normal equations: 

TZro' = n-A + Efc^ (4) 

o=l o=l 



and iExoVo' = aZZxq + bZ!(x<>) 2 (5) 

for B by Gaussian elimination and for A by substitution in (A). Exponen- 
tiation of A and B gives a and b for output. 



LEAST SQUARES FIT: (POWER FUNCTION) 

y = axk-^-ln y = In a + b(ln x). Let In y = y', In a = A, In x = x'. Hence 
y' = A + bx*. To fit a line of this form to a set of observed data (xo,yo), 
we perform the transformations xo' = In xo and yo 1 = In yo. Then we solve 
simultaneously for A and b the nontal equations: 



n 

y ' = n-A + bV"xo' 
^1 

^Txo'yo' = aZLco' + b^(xo«) 2 



(6) 



and Z_Xo'y ' = AZ_Xo« + bZ-kb'r (7) 

for b by Gaussian elimination and for A by substitution in (6). Exponen- 
tiation of A gives a for output. 



COMPUTATION Ftp. STANDARD ERROR, PLOT BACK AND FUNCTION EVALUATION: 

Once parameters have been computed, machine can compute directly y c , the 
plot back value. The deviation, yd = y c - y c , is readily available for 
typeout and for squaring and suraning in computing the standard error. 
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Interpolation and function evaluation are also done directly bv FORTRAN 
using dummy y Q values. Dunmies do not appear as output. 

ERROR ANALYSIS: 

Interrelatedness of summations used for least souares fit makes general 
statements about error analysis difficult. Individual problem approach 
seems best regarding reliability of comnuted parameters. 



PROGRAM CREDITS: 

Southern Services, Inc., Birmingham, Alabama, is contributor of this program. 

Inquiries and comments should be addressed to E. J. Orth, Jr., director, 
Computing Center, at the above address or telephone: Birmingham FA3-5341, 
Ext. 2339. He contributed much to this program. 

We wish to thank Dr. R. E. Wheeler, head of the mathematics department, 
Howard College, Birmingham, for his constructive comments and Mr. Jim 
Moore, IBM, Birmingham, for permission to include his FORTRAN Compressor. 

We wish also to thank our programming colleagues, Mrs. Mary T. Shannon and 
John L. Redmond, for their assistance. The final manuscript is the careful 
work of typists, Mrs. Tallulah S. Armistead and Miss Betty Ross Armstrong. 
Our machine room supervisor, Jim Ewing, says he understands the Operating 
Instructions, so hopefully we submit the program. 

Wade A. Norton, author 
Birmingham, Alabama 
25 July 1962 
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INPUT 



CARD INPUT: 

All card input is one point (observation) per card in the form: 

where I is fixed point identification number processed by 1 620 and must 
appear entirely in columns numbered < 72. Form of Xq and y is floating 
point, End-of-file card contains dunmy (xo»yo) and negative value for I. 
Dunmy y values are used for function evaluation and interpolation, thus 
making all card input compatible. 



TYPED INPUT: 

Under Option 5, two (or three) floating point values are keyed-in at con- 
sole typewriter: 

a RS 
b RS 

c RS (c_ not accepted unless hyperbola). 
OUTPUT (ALL TYPED) 



FUNCTION CODE: Fixed point digit 1 , 2, or 3 on line by itself indicates 
subprogram in use: Hyperbolic, exponential, or power, respectively. 

CURVE FIT (Pass 1 ) : 

Line 1:1, a, b, c (5 only for hyperbola) 

Line 2: a, (In b) (Line 2 only for exponential) 

where I is absolute value of fixed point identification number appearing 
in end-of-file card, and a and b (and if c_) are computed parameters for 
curve. For exponential fit, value a also represents the y-intercept (for 
x - 0) of the linear (semi-log) form, and (In b) the semi-log slope. 

PLOT BACK (Pass 2): 

*o» v o» y c> y d 

where Xp and 2o are tne observed values, yjj is the computed y_ for fitted 
curve, and yd = y c -y c . 

STANDARD ERROR (Pass 2): _ T - 

n, SE, I, ^-(y d ) 2 

where SE is floating point computation of standard error and others pre- 
viously defined. 

INTERPOLATION (-Pass 3) AND FUNCTION EVALUATION: 

where variables indicated are same as described above. "Observed 1 ' y is 



S 



not printed, as it is a dummy used only to make all input compatible. 

PROGRAMMED ERROR MESSAGE (Pass 2, identified bv all fixed point output with 
last negative): 

N, M, -I 

where N is first pass count, M is second pass count, and I is negative, 
.just as read from end-of-file card. A meaningful standard error can be 
computed only on data used to compute parameters. Program utilizes 
three-pass system: (1) compute parameters, (2) compute Standard error and 
plot back, and (3) interpolate. By using three passes, we place no inhe- 
rent limitations on n, number points used in fit. The 1620 cannot be 
sure that same data are processed on Passes 1 and 2 but can be sure that 
same number of cards N, M are processed during respective passes before 
computing standard error. When same number of points are not processed 
each pass, standard error is not computed, and first and second pass counts 
respectively are typed. In one sense this is NOT NECESSARILY an error 
message, as user may not intend to use same number of points each pass. 
(See Option 4, page 16). Number of interpolated values is independent 
of fit, so no count of Pass 3 points is kept. 



FORTRAN ERROR MESSAGES: On Pass 2, a zero divisor will result in FORTRAN 
message ERROR E7 while computing the standard error when n = 3 for the 
hyperbola or n = 2 for the exponential or power functions. This does not 
affect the plot back, but the standard error itself should be ignored. 
The program is subject to other FORTRAN error messages, but none should 
occur in normal operation. Should the messages ERROR Ex or ERROR Fx 
(where x is a digit) appear on the typewriter, consult the FORTRAN write- 
up. Use of fewer points than specified above results in insufficient 
data to define uniquely the curve in question, and a variety of error 
messages of this sort are likely. Under such conditions the output 
should be ignored, but the program remains intact. 
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SAMPLE PROBLEMS 



SAMPLE PROBLEM #1: HYPERBOLA 

The 1620 was given the 12 points listed in the plot back for this nroblem. 
These points are known to lie on both continuous curves comprising the hy- 
perbola (x+7)(v+2) = 30. Results are obvious from glance at plot back and 
interpolation. 

SAMPLE PROBLEM #2: HYPERBOLA 

The 1620 was given the 8 points listed last in the plot back for SAMPLE 
PROBLEM #3. These points are known to lie on both continuous curves com- 
prising the hyperbola xy - 12. Neither plot back option nor interpolation 
option was used. Success is obvious. 

SAMPLE PROBLEM #3: HYPERBOLA 

ThV 1620 was given the 20 points comprising the data for SAMPLE PROBLEM #1 
and SAMPLE PROBLEM #2. Computed values were checked by desk calculator. 
Standard error is high, because points do not cluster about a single hyper- 
bola. 

SAMPLE PROBLEM #4: HYPERBOLA 

The 1620 was given the parameters a = 2, b * -3, c = 360. Function eval- 
uation for Xj, - 4, 6 gave y c = 177, 87 respectively. 

SAMPLE PROBLEM #5: EXPONENTIAL 

The 1620 was given the 22 points listed in the plot back of this problem. 
These are really two sets of points: eleven points on V£ " 2e x and eleven 
points on yo » (e*)/2. Observation of output shows program fitted yi = e x 
to accuracy of input. Standard error is high, because test data fail to 
cluster about single curve. 

SAMPLE PROBLEM #6: EXPONENTIAL 

The 1620 was asked to evaluate y = 3e x for Xq = 0, 1, 2, . . ., 10; ie., 
a - 3 and b - 2.7183 were keyed-in at console typewriter. Results compare 
favorably with desk methods based on identical "input". 

SAMPLE PROBLEM #7: POWER FUNCTION 

The 1620 was given the 10 points listed in the plot back for this problem. 
Five of these points lie on ex 2 and five on exjL Observation of output 
shows pro gram fitted ex*-**? to accuracy of input. Standard error is high, 
because data fail to cluster about a single curve. 

SAMPLE PROBLEM #8: POWER FUNCTION 

The 1620 was asked to evaluate y - 2x? for Xq values listed in typeout; ie., 
a » 2, b - 3. Output is slightly off because exponent is floating point, 
necessitating use of subroutine by FORTRAN. 



c 
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SAMPLE PROBLEMS 



LISTING 



Sample Problem #1, CS 



10 



I0~~ 



-7.QQQQQQQ 



-2.0000000 



30.000000 



-5 . OQQQQQQ 



13.QQQQQQ 



13.000000 



■00000000 



-2.0000000 
-1.0000000 



4.0000000 

3.0000000 



4.0000000 
3.0000000- 



.00000000 
.00000000 



3.0000000 
8.0000000 



1.0000000 
.OOQQQQQO 



1.0000000 
.00000000 



.00000000 

,00000000 



23.000000 
-9.Q0QP0QQ 



-1.0000000 
-17.QOOOOO 



-1.0000000 
■I7.nnnnnn 



.00000000 
.00000000 



-12.000000 
-13.000000, 



-8.0000000 
-7.0000000 



-8.0000000 
-7.0000000 



.00000000 
.00000000 



-17.000000 

-37.000000 



-5.0000000 
-3.0000000 



-5.0000000 
-3.0000000 



.00000000 
.00000000 



53.000000 
12 



-1.5000000 

.00000000 



-1 .5000000 



.00000000 
.00000000 



4.0000000 
6.0000000 



^72727270 
.30769230 



1 



Sample Problem #2, Csj j' 



3333 -.00000000 .00000000 12.000000 



.00000000 3333 .00000000 



1 
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4910312 


Manual Branch to Exponential - 


12 
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Sample Problem #5 


, cs . 10 












1111 


1 nnnnmP. 


■1 — - — -19 

L . / IPZ/jU 












1.0000108 


.99999752 








Sample Problem #3 , CS °P, 




.00000000 
1.0000000 


.50000000 
1 .3592000 


1.0000108 
2.718*04? 


-.50001080 

-1,3591043 










2.0000000 
3.0000000 


14.778200 
10.043000 


7.3890988 
20.085601 


7.3891020 
-10.042601 


1 








4.0000000 
5.0000000 


109.19600 
74.200000 


54.598188 
148.41288 

— * ~ w * ^ B WW 


54.597820 
-74. 2 12flflQ 


8888 


_c 99091 f A 


-1 869 3ii<;4 


37.272s*;* 


6.0000000 
7.0000000 


806.86000 
548.30000 


403.42705 
1096 62*56 


403.43295 
-545.37560 


-5.0000000 


1 % oononn 


167 4P.P.C6 


-15^.38856 


8.0000000 
9.0000000 


5962.0000 
4051 6000 


2980.9299 
8102 987* 

MlUtt JU/J 


2981.0701 
-4051.3873 


-2.0000000 
-1.0000000 


4.0000000 
3.0000000 


9.7052210 
6.9625702 


-5.7052210 
-*. 9625 702 


10.000000 
.00000000 


44052.000 

2 0000000 


22026.148 
1 OOOOlOfl 

— I m 1/ WWW 1 wO 


22025.852 
. 99998920 


3.0000000 
8.0000000 


1.0000000 
.00000000 


2.6649117 
.95001600 


-1.6649117 

-.95001600 


1.0000000 
2.0000000 


5.4366000 

3 6946000 

3 - " J— tjumf 


2.7183043 

7 *8909flft 


2.7182957 

_3.69449P.fl 


23.000000 
-9.0000000 


-1.0000000 
-17.000000 


-.54857030 
-11.730361 


-.45142970 

-5.2696*90 


3.0000000 

4.0000000 


40.172000 
27 299000 


20.085601 
54 598188 


20.086399 
-27.299188 


-12.000000 
-13.000000 


-8.0000000 
-7.0000000 


-7.3669431 
-6.6602919 


— * ♦■— E» W,^ W 

-.63305690 
-,^3970810 


5.0000000 

6.0000000 


296.82000 
201 .72000 

*- v • * i *-WWW : 


148.41288 
40* 42705 


148.40712 
-201.70705 


-17.000000 
-37.000000 


-5.0000000 
-3.0000000 


-5.0334558 
-3. 04218*6 


3.3455800E-02 
4.2l8*600F-02 


7.0000000 

8.0000000 


2193.2000 
1490.5000 


1096.6256 

2980.9299 


1096.5744 
-1490.4999 


53.000000 
4.0000000 


-1 .5000000 
3.0000000 


-1.2291459 
2.1731380 


-.27085410 
.82686200 


9.0000000 
10.000000 


16206.200 
11013.000 


8102.9873 
22026.148 


8103.2130 
-11013.148 


6.0000000 
3.0000000 


2.0000000 
4.0000000 


1 .4525665 
2.6649117 


.54743350 
1.3350883 


22 


5921.7808 


1111 


7.0134993E+08 


2.0000000 
12.000000 


6.0000000 
1.0000000 


3.2929067 
.29512010 


2.7070933 
.70487990 


1 .5000000 

9.5000000 


4.4817204 
13359.551 






1.0000000 
24.000000 


12.000000 
.50000000 


4.1228226 
-.59377110 


7.8771780 
1.0937711 


2 









-4.0000000 -3.0000000 28.676629 -31.676629 
20 38.343302 8888 24993.553 



Sample Problem #4. CS1 



01 
0& 



±2«_L 
-3.K 
36QJS 



4.0000000 
6.0000000 



177.00000 
-87.000000 



3JK 



2.7183K 

.00000000 



1.0000000 
2.0000000 



3.0000000 
4.0000000 



5.0000000 

6.0000000 



7.0000000 

8.0000000 



Sample Problem 



CS 



01 
-00- 



3.0000000 



8.1548991 
22.167460 



60.257802 
163.79877 



445.25415 
1210.3343 



9.0000000 
10.000000 



3290.0514 
8943.3459 



24310.695 
66 083.757 



4912$96 Manual Branch to Power Function 
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S&aqple Problem #7 






2222 


2.7182994 


1.2500001 




1.0000000 


2.7183000 


2.7182994 


6.0000000E-07 


if. ooooooo 

9.0000000 


5 .4366000 
220.18230 


15.377024 
42.374096 


-9.9404240 
I77.flnfi9f 


lo. 000000 
25.000000 


10.673200 
1698.9375 


86.985603 
151.95760 


-76.112403 
1546.9799 


1.0000000 
4.0000000 


2.7183000 
43.492800 


2.7182994 
15.377024 


6.0000000E-07 
- — 28.115776 


9.0000000 
16.000000 


8.1549000 
695.88480 


42.374096 
86.985603 


-34.219196 
608.89920 


25.000000 
10 


13.591500 
593.98197 


151.95760 
2222 


-138.36610 
7822519. 


1.5000000 


4.5124393 






9.5000000 


45.336904 
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<u»pi« P roblem #8. CS ^ \ 



2*B 

.00000000 
1.0000000 
1.5000000 



.00000000 



2.0000000 
6,7499996 



2.0000000 
f.QPOOQQQ 



15.999999 



4.0000000 
5.0000000 



127.99997 
249.99998 



C 



COMPILATION INSTRUCTIONS 



Since the program is supplied in source form, it will be necessary for the 
user to compile the Drogram himself. Note that input/output is designed 
for pre-Format FORTRAN. The source deck i3 numbered serially in columns 
78-80. To assist in the location of cards in the source deck, the first, 
last, and every tenth statement in the program listing have been numbered 
at the right margin. This number corresponds to the serial number of the 
card in which the statement appears. 
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OPERATING INSTRUCTIONS 



Clear Core. 



Set Check Switches: Overflow to Program, 
I/O to Stop, 
Parity to Stop. 

Set typewriter tab stops 17 spaces apart for FORTRAN output. 
Load Program. 



Set CS for desired option as explained below. 

Branch to subprogram: 07500 for hyperbola,* 
10312 for exponential, 
12596 for power. 

The appropriate function code will be printed. 

When the Reader-No-Feed light comes on, read in data with end-of-file card. 
Or, according to option chosen, if typewriter enables, enter function 
parameters (see page 7). 

After subprogram is executed, another subprogram may be branched to. 

Console Switch Settings** 

Option 1: Compute Equation Parameters Only: 0001 

Option 1 allows batch processing. Each set of observations 
nist have an end-of-file card; and after each set is nrocessed 
the function code is printed. When desired set(s) of data has/ 
have been processed, the user raay reset the console switches 
(if desired) and branch to another subprogram. 

Option 2: Compute Eouation Parameters, Plot Back, Standard Error, and 
Interpolation: 1010 

Only one data file mav be processed at a time, since two passes 
through the Reader must be made. On first pass, one or two lines 
of curve fit output will appear. A Reader-No-Feed light also 
signals the beginning of Pass 2, where plot back and standard 
error are computed. After data file is processed by Pass 2, in- 
terpolation points should be entered in the Reader. Interpolation 
points consist of Xq, dummy y Q , and I. Dummy y values and I are 
not printed. Interpolation cards must have their own end-of-file 
card. After interpolation data are read, the function code is 
printed. This is the time to branch to another subprogram (after 
setting switches as desired). 

If the card counts for Passes 1 and 2 are not equal, an error 
message will print after the end of Pass 2 in the form of 3 num- 
bers: number of cards in Pass 1, number of cards in Pass 2, and 
-I. To correct, see note 3, Console Program Switches, page 17. 

* See Comments, page 18. 
** See Note 9, page 26. 
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Option 3: All Features of Option 2 except Interpolation: 0010 

After Pass 2 the function code is printed. This is time to 
branch to another subprogram (after setting switches as desired). 

Option 4: Compute Parameters, Some Plot Back, Some Interpolation: 0010 

Option 4 is a variation on Option 3 and uses the same switch 
settings. Procedure is the same except as noted below. Option 
4 will be applied under these conditions: 

(a) standard error and sum of squared deviations are not 
desired, 

(b) only partial plot back desired, 

(c) interpolation is or is not desired. 

The data file for Pass 2 consists of all points to be plotted 
whether plot back or interpolation. Pass 2 data file must have an 
end-of-file card. On interpolation, y Q and y^ of plot back will 
appear and should be disregarded. Also ignore the last line of 
print before the typing of the function code: it will be either 
an erroneous standard error or an incorrect error message. Option 
4 ends after typing of the function code. When the halt (48 in op 
register) occurs because of false error message, turning CS3 OFF 
and hitting start causes tvping of the function code. 

Option 5: Function Evaluation: 0100 

This option is especially valuable when interpolation or plot back 
is "after-thought''. When 37 op code appears in console op register, 
key-in function parameters. When Reader-No-Feed light appears, load 
data file which consists of data cards (containing xo, dummy y , I) 
and end-of-file card (where both Xq and y are dunray values and I is 
negative). After the evaluation is completed, the function code is 
typed. 

CONSOLE PROGRAM SWITCHES (DETAILED DISCUSSION) 
All comments below refer to use of a single subprogram. 

1. CS2. In Pass 1 of Options 1 through 4, when Reader-No-Feed light 
appears on console, CS2 has been interrogated. It may be turned ON 
now or at any time before last pass is completed, if next run is to 
be function evaluation. When using Option 5, best olace to turn CS2 
OFF (if next run is fit) is at Reader-No-Feed before type-out. After 
completing any option, machine tynes function code and goes to one of 
two program points. With CS2 ON, 37 appears in op register to indicate 
function evaluation block has been entered. With CS2 OFF Reader-No- 
Feed light appears to indicate Pass 1 block has been entered. 

2. CS3. Plot back can be turned OFF or ON while in progress. An 
"Option 2a" involving three passes will vield parameters, limited 
plot back, standard error, and interpolation. Settings and steps 
are as for Option 2 with following modification: (a) Arrange observed 
data file, so points for which plot back desired are at front of file 
on Pass 2, (b) Turn CS3 OFF while last point of desired plot back is 
being typed, (c) Use Pass 3 for interpolation. Such "Option 2a" 
could replace Option 4. Op 4 has advantages of less machine time and 
requiring no change of CS settings during run. Op 2a has advantages 
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of consistency with Op 2, valid standard error furnished, and no 
dummy y typeout for interpolation. 

CS3. After "ERROR MESSAGE" (N, M, -I) machine halts. If true 
error has occurred, turning CS3 ON permits repetition of Pass 2: To 
repeat Pass 2, clear Reader, locate error, load Reader-Hopper with 
correct data file and hit Start. Remember CS3 controls plot back, 
so when repeating Pass 2, change CS3 setting, if necessary, before 
hitting Reader-Start. With CS3 OFF program returns to type function 
code and interrogate CS2. 

CS4. CS4 is interrogated at the end of Pass 1. If CS4 is ON, CS2 
is interrogated as stated a bore. If CS4 is OFF, the program con- 
tinues to Pass 2. 

CS1. CS1 is interrogated at the end of Pass 2 when no error con- 
dition detected. If CS1 is ON, the program continues to Pass 3 for 
interpolation. If CS1 is OFF, CS2 is interrogated as stated above. 
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C THIS PROGRAM FITS BY METHOD OF LEAST SQUARES CURVES OF FORM 1 

C Y=B+C/(X-A) (AN HYPERBOLA), Y=A(B**X) (AN EXPONENTIAL) , AND 

C Y=A(X**B). LEAST SQUARES METHOD IS APPLIED DIRECTLY IN CASE OF 

C HYPERBOLA AND TRANSFORMATION TO LINEAR FORM (SEMI -LOG OR LOG-LOG 

C RESPECTIVELY) IS MADE IN OTHER TWO CASES. MANUAL BRANCHING VIA 

C 1620 CONSOLE IS NECESSARY TO MOVE FROM ONE TYPE FIT TO ANOTHER. 

C (FOR BRANCHING DETAILS, SEE OPERATING INSTRUCTIONS.) FIRST PART OF 

C PROGRAM IS HYPERBOLIC FIT WHICH FOLLOWS. BRANCHING IS NOT NECES- 

C SARY TO ENTER IT INITIALLY, BUT IS OTHERWISE. 

DIMENSI0NS(2,if) 10 
03 D05K=1,2 

05 S(K.J)=0. 
PRINT, 
PRINT, 1 
PRINT, 
MM=3 

IF(SENSE SWITCH2)1,12 
1 ACCEPT. S(2, k), B.C 

L=2 20 
GO TO 6 
07 IF(I)3.4,4 
k S(1,2)=B+(C/(X-S(2,4))) 
GO TO 1429 
12 L=1 

66 READ.X.Y.I 

GO TO (2.7.51, 7, 1007, 1032,1129, 1129, 2007,2032, 2129, 2129), L 
02 I F( 1)10,8,8 
_08 S(1,3)=S(1,3)+1. 

S(1,1)=S(1,1)+X 30 

S(2,lj-S(2,1)+X**2 

SU,2>S(!,2)+Y 

S(2,2)=S(2,2)+X*Y 

S(2,3)=S(2,3)+Y*Y 

S(2,1»)=S(2,4)+X*X*Y 

S(1,4)=S(1,4)+X*Y*Y 

GO TO 6 

10 B22=S(1 ,2)*S(2,2)-S(1,1)*S(2,3) 
B21=S(2,1)*S(1,2)-S(1,1)*S(2,2) 

B23=S(1.2)*S(2,4)-S(1.1)*S(1.4) 40 
B=B22*(S(2,2)*S(1,2)-S(1,3)*S(1,4))-B23*(S(1 ,2)**2-S(1 ,3)*S(2,3)) 
B=B/(B22*(S(1,1)*S(1,2)-S(1,3)*S(2,2))-B21*(S(1,2)**2-S(1,3)*S(2,3))) 
S(2,4.)-(B23-B21*B)/B22 
"CiT5(Tpi7-S(2,2)*B-S(2,3)*S(2,4))/S(1,2) 

C=C+S(2,4)*B 

l=-l 

PRINT,I,S(2,4),B,C 

PRINT, 

N-S(1,3) 

— TF( SENSE SWITCH^ 3 ,3V 50 
5k L=3 

GO TO 1131 
51 IF( 1)1*9,50,50 
50 M=M+1 

S(1,2)=B+(C/(X-S(2,4))) 

S(I.4)=Y-S(1 .2) 

S(2,1)»S(2,1)+S(1,4)**2 

IF(SENSE SWITCH3)45,6 
45 PRINT,X,Y,S(1,2),S(1,4) 

GO TO 6 60 
49 IF(M-N)52.53.52 
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53 S(2,2)=SQR(S(2,1)/(S(l,3)-3.)) 

GO TO , 1 

60 IF(SENSE SWITCH1)*6,3 L , ,. . 
52 PRINT,N,M,I 

PAUSE ^7 , 

61 IF(SENSE SWITCHf)5*,3 
*6 L«* 

C HERE°BEGINS PROGRAM FOR EXPONENTIAL (SEMI-LOG LINE) FIT, Y-A*(B**X) 70 

1001 D01002K-1.2 
001002J-1 ,* 

1002 S(K,J)-0. 
PRINT, 
PRINT, 2 
PRINT, 

IF(SENSE SWITCH2) 1029, 1006 

1006 L-5 . 80 
GO TO 6 

1007 IF(I)1018, 1008, 1008 

1008 S(1,3)-S(1,3)+1. 

9999 Y-LOG(Y) ■ ' 

S(1,2)«S(1,2)+X 
S(2,2)-S(2,2j+X**2 
S(i.*)«S(1,*)+Y 
S(2,*)«S(2,*)+X*Y 
GO TO 6 

s?2;i>!s!iH)^ 90 

' S( !lRs(i:^-SC2.li*S(1,2j)/S(1.3) 

^h.M-EXPts{l,lD 

I — I . 

" ' ' GO tO (1019, 2019), MM 

1019 Si2 J 3)-EXP(S(2jj) 
PRINT.I.S(2.*).S(2.3) 

PRINT, S (2jj»hsl2»U _.. 
PRINT 

I F( SENSE SWITCH^IIOOIJUO 

1130 L^ 100 

1131 SdJj-O. 
M-0 

GO TO 6 
-1072 17(1718*3,1036,1036 

1036 M-M+1 , , % _. 

5(1,2}-S(2.*)*(S(2,3)**X) 

S(1.*)«Y-S(1,2) 
th*U-S(2,f1 + S<1,^)**2 

I F( SENSE SW I T CH3) 1039^6 
1039 PRIN T , X.V, 5 ( 1 7277507*1 110 

GO TO 6 

10*3 5(2,2)^(5(2, lJ/(S(1,3)-2.)) 
10S5 I--I , , % 

PRINT. N, S(2.2) , I , S(2j 11 

PR | NT 

GO TO* (10**, 20*4.60), MM 
10** 1F(SENSE SWITCH!) 1052, 1001 

10*7 PRINTS, M, I 12Q 
PAUSE , 4 

GO TO (10*8,20*8.61 ),MM 
"T6*ff TF( SENSE SWITCH*) 11 30, 1001 

1029 L-7 /3N 

1030 ACCEPT, S(2,*),S(2,3) 
GO TO 6 

H29 IF(j)1001, 1329.1329 

1329 S(1.2)»S(2.*)* S 2.3)**X) 



1^29 PRINT, X,S(1 ,2) 

GO TO 6 

1052 L=8 20 a WO 

- HERE°BEGINS PROGRAM TO FIT LINEAR LOG-LOG PLOT FOR CURVE Y=A*(X**B) ' " 

2001 D02002K=1,2 
D02002J=1 ,* 

2002 S(K,J)=0. 
PRINT, 
PRINT, 3 
PRINT, 
MM=2 

IF(SENSE SWITCH2)2029,2006 140 

2006 L=9 

GO TO 6 

2007 IF( 1)1018.2008,2008 

2008 S(1,3)-S(1,3)+1. 
X=LOG(X) 

GO TO 9999 
2019 S(2,3)=S(2,1) 

PRINT, l,S(2 ( *) t S(2,1) 
PRINT, 

IF(SENSE SWITCH*)2001 ,2130 150 
2130 L=10 

GO TO 1131 
2032 IF( 1)18*3,2035,2035 
2035 M=M+1 , 

S(1,2)=S(2 | *)*(X**S(2,3)) 

S(1,*)=Y-S(1.2) 

S(2,1)=S(2,1)+S(1.*)**2 

IF(SENSE SWITCH3)1039,6 
20** I F( SENSE SWITCH1)2052.,2001 

20*8 iF( SENSE SWITCH|n 1131, 2001 160 
2029 L=11 S3\ 

GO TO 1030 
2129 IF( 1)2001. 2329.2329 
2329^(1,2)=S(2»*IX**S(2,3)) 

GO TO 1*29 
2052 L-12 

GO TO 6 

END 166 

END OF COMPILATION 
LOAD SUBROUTINE DECK 
THEN PUSH START 



SBCOL TABLE PRIOR TO SPREAD (INCOMPLETE) 
362000000 10048ftS$ This listing was obtained by using the 1-per-card 
38 1 8860 program in the following manner: The last card in 
the panch stacker at time subroutines called for was 

OOOOOOOOOOOOOOOOOOOOgOOOOOgQ12T35320T329gOOOOOgOlT T 

358807052T32 1 6oT035PO0O0OgQloT3 1 28571 3OT2956o7O08gQ 

0000g009'T2868oT006T3 i »ifi»g2'029T2620g2'002T 2 g96g2'00tg00 

OOOgOOB72472gT329T2384gT030gOOOOOg007T3424g?048T234 

0070487256407052734040704472 1 80gT044gl 200000007l 932 

gTo43T22OOo7o47Tl796o7039Tt608o7036Tl864o7843g00000 

gQOOPOOOOOgOQgT l472ffTl 30T3024oT019Tl33207019T0708g9 



999T0672gTO08T09O8gT0 1 8gOOOOOg00570584o7006T2360070 
29T0336gT0O2To^12CTooig0000Qg00yT0260gg06lT0280gg04 
6T01 8QggQ60Tl992gT045gl 3000000gT01 1 2ggQ53T0200gg052 
g9976gg045g400000000g9776gg050T0044ggQ49Tl 496gTl 3 1 g 

9688gg054g5000OOOOOgOOOO00gO6g0O0OO0g03g0OO0OOgO5gQ 
00000g04g2 72 730000^2 72 7 1 OOOOgOOOOOOOO 1 $2 72 72 OOOOgOO 
OOOOObOgl 1 000000008 1 96CgQ08g85 76PC0 1 0T347607l 29T3 1 6 
007032T29000?007T241 6071 29T1 5 5207032T061 6gT007g9720 
ggQ5 1 g8 1 4007k)02g800000000T2 5200T429O*0O0O00gQ2S7O00O 
000Qg7896gOO0M?90OOOCK)00g7840ggOO7g7988g0O06g00OO0g 
00ZT30O0OOOOOg30O0OOO00^2000O0O0Og7964gQ0 1 2g7 772ggO 
0lg0OOOQgOO3T45 4000 00000 oooogooTDoooooooooTi 0000000 
0g20OO000OOg7524o7fo05g7500gg003ggO02T9809g2OO000000 
Tjj 63554600^1 63550OOOg3566246OQg35662O0O0S2495S4600g 
249S50000S258594600g258590000g3564746005356470QOOl?5 



67574600^567570000 

noted, uoapixation was then completed. One-per-card program 
wag then loaded a* far as card noted above. Appropriate place for 
type-out was determined, and this listing was produced. Purpose 
was to determine beginning of exponential and power subprograms 
which are underlined. 



o 



© 



22 

MAINTENANCE 



LIKELY CHANGES 

Programmed CS control should give sufficient versatility to meet most 
needs for these curves, especially at installations where machine con- 
figuration is compatible with that for which this program was designed 
Source program cards numbered serially in columns 78-80 are furnished. 
Most likely changes will be to make I/O compatible with user's needs. 

However, should alteration after compilation prove desirable, the last 
page of PROGRAM LISTING is tvpeout of most of svmbol table prior to 
spread. Details are contained on that page. 

In this regard, compressed program (p. 26) is more easilv altered, as 
an absolute address appears in columns 1-5 of each card. 
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1. The symbol ZZ in the writeuo is equivalent to ZZ in the FORTRAN listing 
and was so written to eliminate any confusion arising from ZZ incorrectly 
implying a oroduct to some readers. 

2. When parameters or variables, a, b, etc., appear in text, they are under- 
lined for reading ease, When they appear in equations or. parentheses, 
they are not. 

3. RS means perform release and start operations at 1620 console. 

L. AK INTERPRETATION OF THE STANDARD ERROR 0TH2R THAK TESTING THE SIQCIF- 
IGAKCF CF AK ADDITIONAL PARAMETER: This discussion also pertains to poly- 
nomial curve fitting, where use of the standard error in testing additional 
parameters is common. See "rogram Number 7.0.002 (Graves). 

The given empirical data quite likely do not lie on any single plane 
curve, because they are not the function of two but several variables. 

However, if the data represent the function of one main independent vari- 
able with the variability in the dependent variable due to lesser impor- 
tant 5 '" i"~f"-"''>nt variables, the curve fitting may possibly be treated 
statistica_Uy . 

An extremely useful assumption (provided it is valid, of course) is that 
the given data represent a randan sample of y's from a normal population 
Y which has constant sigma square about the computed value y^. given by 
y c = f(x) where f(x) is the fitted function. It is useful because of 
its ease of computation and ease of interpretation. 

The assumption of constant sigma square for Y irrespective of x is in par- 
ticular the assumption that errors of measurement and the effects of other 
independent variables of considerably less importance than the main one 
yield an aggregate that is normally distributed with constant sigma square, 
and independent of x. 

Describing this latter variability then is the real issue. 

What we would like to say is that approximately 95$ of the population of 
values which the true function of several variables will take on lie be- 
tween y c -2sigma and y + 2sigma for all x. That is instead of fitting 
only the single curve y c , we would fit in reality the family of curves 
y c + k(sigma), which have statistical significance and for which the 
standard error is a valuable estimate of sigma. 

Our question then is, "Is the variability of the less important variables 
normally distributed with constant sigma square?'' If so, we can make a 
reasonable assumption that it is true for y in Y. Although the validity 
of such an assumption is beyond the scope of this writeuD, it should be 
based uoon more than the desire of the user for "nice" results. 

When statistical reliability is not needed, the sum of the deviations 
squared can be studied directly. This value is also furnished. 
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5. '.ftien fitting a curve to points, order in which points are given to 1620 
does not matter. Order on Pass 1 need not be same as for Pass 2. How- 
ever, if standard error desired, Pass 2 must use all points. 

6. Numerical analysts are in general agreement that when a curve is fitted 
to emniric?l data, its use for extrapolation is tc be approached with 
great caution. Interpolation is much more reliable. 

7. This program may be compressed to 203 cards by use of General Purpose 
Dump and FORTRAN Compressor written by Jim Moore, IBM, Birmingham. Use 
of the Comnressor (Anpendix C) requires Indirect Addressing. 

3, If manual branching from subprogram tc subprogr.'ir.i occurs only at points 
indicate-i in operating instructions, result is to "frame" any work done 
by given subprogram between tynings of its subprogram code number. 

9. With reference to console switch settings: The notation XXXX (X = 0, 1) 
refers to CS1 , CSS, CS3, CS4 from high to low order digit, respectively; 
means OFF, 1 ON. 

10. Similar treatment could be given to the general logarithmic equation 

(y = b(ln x)+ In a), several other forms of hyperbola, and several forms 
for the circle and ellipse. Another possibility (requiring solution of 
a 6x6 system) is to fit a curve of the form ax 2 + bxy + cy z + mx + ny + 
k = 0, the general second degree curve. A 40K program similar to this 
20K program could include this one and contain several other such sub- 
programs . Likely the general equation would require a separate program, 
but still other subprograms from the list above could be packed into a 
second 40K program. Each curve type in this 20K program and in the above 
list has the advantage of tending not to be too sophisticated for ordi- 
nary empirical data to 3 or 4 significant figures. Southern Services 
plans to develop programs for these other curves as the need arises. 

1 1 .This program leaves unutilized approximately 5200 of the basic 20,000 
core locations. This should be more than sufficient to permit addition 
to this 20K program of one more three-parameter subprogram of the user's 
choice. 
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COMPRESSES FORTRAN DECKS TO LESS THAN HALF ORIGINAL SIZE. EVEN COMPRESSED 
SPS DECKS MAY BENEFIT. ALLOWS SNAPSHOTS OF MEMORY OVER ANY INTERVAL. ROUTINE 
AND DUMPED DATA ARE SELF-LOADING WITH AUTOMATIC BRANCH TO PROGRAM. TABLES AND 
PROGRAM OVER POSITION 100 ARE INTACT IMMEDIATELY AFTER DUMPING AND ALSO UPON 
RELOADING. REQUIRES IDA. 

INSURE COLUMNS 56-60 ON CARD 10 CONTAIN HIGHEST 1620 ADDRESS WITH FLAG. 
LOAD. INSERT DUMP STARTING ADDRESS WITH FLAG. R-S. INSERT LAST ADDRESS TO BE 
INCLUDED WITH FLAG. R-S. INSERT PROGRAM STARTING ADDRESS. R-S. INSERTS ARE •£ 
OPTIONAL WITH ifOO, X9999, AND 7500 ASSUMED AS LIMITS. IN ANY EVENT POSITIONS _ 
i»00 AND 401 ARE INCLUDED TO RETAIN RECORD MARKS, IF ANY. COLUMNS 1-5 ON DUMPED ^ 
CARDS DEPICT LOCATION OF DATUM IN COLUMN 6. DUMPS WHETHER OR NOT TABLES ARE IN g 
MEMORY. 1 _ 

FORTRAN USERS SHOULD DUMP TWICE IF THERE IS MUCH ROOM BETWEEN THE PROGRAM <_o 
AND THE SYMBOL TABLES, NOTING THE BOUNDS OF EACH. THE TWO DECKS ARE THEN COM- 53 
BINED REMOVING THE 4 REAR CARDS FROM THE FIRST DECK AND 5 FRONT CARDS FROM THE » 
SECOND DECK. OTHERWISE, DUMP ALL OF MEMORY AND REMOVE CARDS WITH ZEROS FROM oz 
COLUMN 6 TO 80. So 

TO ELIMINATE THE HALT BETWEEN LOADING AND PROGRAM START CHANGE THE HALT 
COMMAND ON THE THIRD CARD FROM THE REAR TO NOP AND THE ADJACENT ALFA DOUBLE 
DIGIT REPRESENTATION FROM GO TO GONE WITH ALFA RECORD MARK. »S 
NUMERIC BLANKS WILL BE DUMPED AS BLANKS. LOCATE SAME AND PUNCH WITH 4-8. 2£ 
36000 00O5O036OO3OOO05O038003O00O4OO360O3O0OO5O0380O3O0OO40O36O0320OO5OO38OO322O o 
O4O036O000O0O5004900O0000O0000000O00O00O0000OO000000O00O00000000000O0000O00O000O 
36OO3O000500360O3220050036000G00050049O00S20O00O0000O0OO00O00O0000000O00OOO0O00O 2> 
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